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CENTRAL CONFIGURATIONS AND TOTAL 
COLLISIONS FOR QUASIHOMOGENEOUS N-BODY 

PROBLEMS 

Abstract. We consider n-body problems given by potentials of 
the form + ^ with a,b,a,/3 constants, < a < b. To analyze 
^vq | the dynamics of the problem, we first prove some properties related 

_ . to central configurations, including a generalization of Moulton's 

qj'^ theorem. Then we obtain several qualitative properties for colli- 

. sion and near-collision orbits in the Manev-type case a — 1. At 

I the end we point out some new relationships between central con- 

. figurations, relative equilibria, and homothetic solutions. 



Oh- 



1. Introduction 



The n-body problem studied here is given by a potential of the form 
4 + 4; where r is the distance between bodies and a, 6, a, f3 are con- 
stants, < a < b (see [U EE])- In the first part of the paper we treat 
the general problem, and in the second part we focus on the case a = 1. 
^ . The function f^ + ^z, called quasihomogeneous because of being the sum 

of homogenous functions of different degrees, generalizes classical po- 
tentials, such as those of Newton, Coulomb, Birkhoff, Manev, Van der 
Waals, Libhoff, Schwarzschild, and Lennard- Jones. Thus, the applica- 
^ . bility of the quasihomogeneous n-body problem ranges from celestial 

Q\ \ mechanics and atomic physics to chemistry and crystallography. 

Although many properties of the Newtonian n-body problem have a 
correspondent in the homogeneous case, this is not true for nonhomo- 
geneous potentials. On one hand, the transposition of known results is 
far from trivial; on the other, new properties show up. 

An intriguing aspect we will point out in this paper refers to central 
configurations, which are crucial for understanding the dynamics of 
the n-body problem (see [13]). The central configurations of the quasi- 
homogeneous potential are in a certain relationship with the central 
configurations of the homogeneous functions that form this potential. 
Thus, we will introduce here the notion of simultaneous central config- 
uration and will investigate its connection with the classical concept. 
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In Section 2, we define the quasihomogeneous n-body problem and 
write down the equations of motion. In Section 3, we introduce the 
concepts of central and simultaneous central configuration, the latter 
being specific to quasihomogeneous potentials. Section 4 deals with 
collinear central configurations. Using critical point theory, we prove 
a generalization of Moulton's theorem by showing that the number of 
collinear central configurations of n bodies is n\/2. Starting with Sec- 
tion 5, we restrict our study to Manev-type problems, [5], i.e. those 
given by potentials of the form - + -4, and show that there are ex- 
actly two planar central configurations in the 3-body case. Section 6 
introduces a framework for the study of collision and near-collision or- 
bits, which is performed in Sections 7 and 8. We study in detail the 
network of collision solutions and determine the relationship between 
central configurations, on one hand, and relative equilibria and homo- 
thetic orbits, on the other hand. It is important to note that if in the 
homogeneous case the correspondence between central configurations 
and homothetic solutions is one-to-one, this fails to be the case in the 
quasihomogeneous problem. The relationship between central config- 
urations and relative equilibria remains unchanged, i.e. one-to-one, in 
the quasihomogeneous case. For Manev-type potentials, homothetic 
orbits are less likely than in the Newtonian case, in the sense that they 
show up only for simultaneous central configurations. 

2. The Quasihomogeneous n-BoDY Problem 

We will start with defining the planar quasihomogeneous n-body 
problem. Consider the linear space 

n 

(1) n = {r = (n, . . . , r n ) G (R 2 ) n \ ^m,r, = 0}, 

8=1 

where m* > 0, i — 1, 2, . . . , n, are the masses of the n bodies and r iy i = 
1,2, ... ,n, represent their coordinates. Notice that Y^i=i m i r i = fixes 
the centre of mass at the origin of the coordinate system. Let 

(2) Aij = {(n, . . . , r n ) G Q\Ti = rj}- A = (J A tj . 

We call A the collision set. The potential U of the system is a function 
defined on the configuration space Q = Q \ A and is given by 



u = w + v, 
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where W is a homogeneous function of degree —a, a > 0, 
(3) W(r h ...,r n ) =^1^ 

and V is a homogeneous function of degree —b, b > a, 

The equations of motion of the n bodies define a vector field X on 
the tangent bundle T(Q). The configuration space of the system is Cl 
and the cotangent bundle is T*(Q). Let p = M _1 r be the linear 
momentum of the system of particles, where M is the diagonal ma- 
trix M = diag (mi, m 1; m 2 , m 2 , . . . , m n , m n ). Then the equations of 
motion can be written as a Hamiltonian system, 

dH 

(5) 1 ^ 

v ; dH 

P = ~ IT - ' 

where if : T*(Cl) — > 1R is the Hamiltonian function given by 

(6) if(r,p) = ^p'M^p - C/(r). 

Here T = | p* M _1 p is the kinetic energy. The total energy if is a 
first integral for the system ([51); this means that T — U = h (constant) 
along any orbit. Other integrals are given by the linear momentum, 
Y^i=i m i^h an d by the angular momentum, J : T — > R, defined as 



(7) J(*,v) = Yl 



a 

' TUiTi X Vj. 

i=l 



Notice that the relationships for the centre of mass, Y^=i m « ri = 0> anc ^ 
linear momentum, 5^t=i m i^i — 0, together with the energy integral, 
T — U = h, reduce the dimension of the Hamiltonian system (jHJ) from 
An to An — 5. We also introduce the scalar product, 

(8) (r, f ) = r*Mf , 

which allows us to write the moment of inertia as 



(9) J=(r,r) = ^ 



n 

2 



i=l 
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3. Central Configurations 

Central configurations play a crucial role for understanding the dy- 
namics of n-body problems [13] . In particular, they have led to impor- 
tant theoretical investigations, such as Saari's conjecture, which has 
remained open for more than three and a half decades [6] , and are con- 
nected to Smale's 6th problem [15], originally proposed by Wintner in 
1941, [16] (see also [91 H]). In this section we will define central con- 
figurations and analyze the particular aspects this concept encounters 
in the quasihomogeneous case. 

Definition 1. A configuration r 6 Q is called central if there is a 
constant a such that 

(10) W(r) = trW(r). 

Using the fact that the functions W and V are homogeneous of 
degree —a and —b, respectively, and applying Euler's theorem for ho- 
mogeneous functions, we find that 

-aW(r) - bV(r) 

ff = — m — ' 

Definition 2. We call r 6 Q a simultaneous central configuration for 
the potentials W and V if there are constants o~\ and 02 such that 

VW(r) = 0iVI(r) and W(r) = a 2 V/(r). 

Using the fact that W and V are homogeneous functions of degree 
—a and —b, respectively, we find that 

/ 10 n -aW{r) -bV(r) 

(12) ox = — —, \ and a 2 — 



2J(r) 2/(r) 

Note that if r is a simultaneous central configuration for W and V, 
then r is also a central configuration for U = V + W. The converse is 
not necessarily true. 
Let 

S Io = {vett\(r,r) = I Q } 

be the sphere relative to the metric given by the scalar product, and 
denote by 

Sl = S Io \A = {ven\(r,r) = I } 

this sphere minus the collision set. Then the central configurations 
with moment of inertia Iq can also be defined as the critical points of 
Usj, where Usj '■ Sj — > M. is the restriction of the potential U to S} . 
Denote by C n the set of central configurations of the quasihomogeneous 
n-body problem. 
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Definition 3. We say that two relative equilibria in Sj are equivalent 
(and belong to the same equivalence class) if they can be made congruent 
by the induced S 1 action on S} , that is, if one is obtained from the other 
by a rotation. 

Let C n denote the set of equivalence classes of central configurations. 
Note that this definition differs from the one used in the Newtonian case 
(see [HE]), where two central configuration are called equivalent when 
one can be obtained from the other by a rotation and/or a homothety. 
This change is necessary in the quasihomogeneous case because the set 
C n is invariant under the action of the group S 1 , but not necessarily 
under the action of homotheties (see Section 7). 

Clearly, I and A are invariant under the action of S . Thus, we 
can conclude that Sj is diffeomorphic to the (2n — 3)-dimensional 
sphere S 2n ~ 3 (which is actually an ellipsoid E 2n ~ 3 ) with all the points 
A removed, that is, 

S* o = E 2n ~ 3 \ (E 2n - ?J n A) « S 2n ~ 3 \ (s 2n - ?J n A). 

Since Ugj is invariant under the action of S 1 , it defines a map Usj '■ 
Sj Q /S l — > R. If we let n n : S} — > S} Q / S 1 denote the canonical projec- 
tion, A = ir(E 2n ~ 3 n A), and recalling that E 2n ~ 3 / S l « S 2n ~ 3 / 'S 1 w 
<CP n ~ 2 (the complex projective space), we are led to investigate the 
critical points of U Sl : CP n " 2 \ A -> R. 

Consequently we can show that the set of equivalence classes of cen- 
tral configurations with fixed moment of inertia Iq is given by the set 
of critical points of the map U $j '■ CP n ~ 2 \ A — > R. More precisely, we 
have proved the following property: 

Proposition 1. For any choice of masses in the planar n-body problem 
with a quasihomogeneous potential, n > 2, the set of equivalence classes 
of central configurations with moment of inertia Jo is diffeomorphic with 
the set of critical points of the map Ug I : CP n ~ 2 \ A — > R. 

4. Moulton's Theorem for Quasihomogeneous Potentials 

We will now study collinear central configurations and, using critical 
point theory, will calculate the number of classes of such configurations 
for any number n of bodies. The goal of this section is to prove the 
following result, which generalizes a theorem obtained by Forest Ray 
Moulton in 1910, HH]. 

Theorem 1. For any choice of masses in the n-body problem with 
a quasihomogeneous potential, U , and any given moment of inertia, 
Iq, there are exactly n!/2 classes of collinear central configurations. 
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In other words, there are n\/2 classes of central configurations r = 
(r x,..., r n ) ; where all belong to the same straight line through the 
origin. 

In preparation for the proof, choose some line / in M 2 . This defines 
a subset Qi C Q of r = (r 1? . . . , r n ) such that each is on the line /. 
Let S t = S Io n Vli and S? = Sj \ n(Sj n A). When S l acts on S Io , only 
the rotation by ix radians leaves S\ invariant. Thus the group Z 2 acts 

-Us 

on Si, and on the quotient we have RP n ~ 2 \ A C CP"" 2 \ A — U R, 
where MP™ -2 = 5 , ;\Z 2 is the real projective space, naturally contained 
in CP n ~ 2 . Here U,s [ is induced by the potential energy. From these 
considerations we obtain: 

Lemma 1. The set of equivalence classes of collinear central configu- 
rations with moment of inertia I is diffeomorphic to the set of critical 
points of U Sl : CP™ -2 \ A -> R that lie in RP n ~ 2 \ A C CP"" 2 \ A. 

So in order to describe the collinear central configurations, it is suf- 
ficient to obtain the critical points of the potential that lie in the real 
projective space. In general, a critical point of a function restricted to 
a submanifold is not necessarily a critical point of the function on the 
ambient manifold. However, we have the following result: 

Proposition 2. If r E RP n ~ 2 \ A is a critical point of U Sl : RP n ~ 2 \ 
A — > R, then r is also a critical point of U Sl ■ CP™" 2 \ A — > R. 

To prove this, we first need to know the derivatives of the potential 
function, which are given below. 

Lemma 2. For given masses mi, ... , m n and U = W + V, 

(1) The first derivative of U : Vt — > R is 

DU(v)(v) 



rriimj 



ri-Tj\ 

rriimj 



\a+2 



K Ti - Tj, Vi - Vj) 



I r • r • 1 1 fr+ 2 



for v G fi. 



The Quasihomogeneous n-Body Problem 



7 



(2) The 2nd derivative is 



^(r)(v,w) = a^ 

a + 2 , 



TfljTflj 



y» 1 1 Qi~ \~ 2 

ii 1 1 7 



-^{Yi - Tj, Vi - Vj)(ri - Tj, Wi - Wj) - (v< - Vj-, Wj - 



+ b > — 

^ r . _ r .116+2 

' \ \W - l W* ^ ~ r i' v * " v i)( r < - r i> w * - w i) " ( Vi " v i' w * 
VII 1 * L j\\ 

where v, w G Q. 
(3) T7ie 2nd derivative of the restriction U : S} — > R is: 

^ 2 C//(^ )(r)(v,w) = D 2 f/(r)(v, w) + flW(r) + 6V(r) (v, w). 

ifere (•, •) denotes the usual inner product in R 2 , || • || the norm in R 2 ; 
and I the moment of inertia. The same formulas are valid in R, R 2 
and R 3 . 

Proof. All the equations above can be derived by differentiating in local 
Cartesian coordinates. □ 

Now we can give a proof of Proposition [2j For Vj G R 2 , let Vj = 
{vj, v^-) where vj G / and vj- G l 1 . Then we can write v = (v", v -1 ) with 
v 11 = (vf, . . . ,Vn), v- 1 = (v±, . . . ,Vn) for each v G Q. If r G Si C S/ , 
r £ A, we have T P (5j ) = {v G fi|(v,r) = 0} and T r (Sj) = {w G 
flj|(w, r) = 0} where, as usual, Q is endowed with the mass scalar 
product. If v G T T (Sj ) and v = (v^v- 1 ), then v" G and (v, r) = 
(vll,r). Thus v" G T r (Sz), because (v, r) = implies (v",r) = 0. 

By Lemma [2] it follows that if r G Si \ A and v G T r {Si ), then 
L>[/(r)(v) = £>C/(r)(v"). So DE/(r)(v") = implies that DU(r)(v) = 
0. This completes the proof of Proposition [2J 

Lemma 3. RP n_2 has n!/2 components. 

Proof. Let r = (ri. . . . , r n ) G Sj \ A and let ri < . . . < r n G R (we 
use the fact that the rj are all distinct). Let a = {i±, . . . ,i n ) be an 
arbitrary permutation of the numbers (1,2, ... ,n). If we apply the 
permutation to the initial vector r, we map it to a different component 
defined uniquely by the given permutation. Therefore the set Si\A has 
n! components and the quotient space RP n ~ 2 \A has 77-1/2 components. 

□ 
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We can now prove Moulton's theorem for quasihomogeneous po- 
tentials. By applying part (2) and (3) of Lemma [21 we see that 
D 2 U / (Si \ A) is a positive definite form, and consequently U is con- 
vex. This shows that U has a unique minimum in each component 
of ILP™ -2 . Thus there are n!/2 critical points and hence n!/2 central 
configurations. 

Remark 1. We have identified the symmetric central configurations, 
otherwise the number of classes of central configurations would be n\. 



5. Planar Central Configurations 

In this and subsequent sections, we will restrict our study to Manev- 
type quasihomogeneous potentials, namely those U for which a = 1 
(see also [5]). They form an important class of quasihomogeneous po- 
tentials, derived from the Manev law, which can explain the perihelion 
advance of the planet Mercury within the framework of classical me- 
chanics (for more details see [I] and [3]). Since for any planar central 
configuration in the Manev-type three body problem, the mutual dis- 
tances are geometrically independent, we can solve the equations defin- 
ing the central configurations in terms of the mutual distances. To be 
precise, we state here a result whose proof can be found in [2]. 

Lemma 4. Let u = /(x) be a function with x = (xi, X2, • • • , x n ), x\ = 

9i(y), x 2 = <fe(y),- ■■, x n = 9n(y), y = (yi, Vi, ■ ■ ■ , y m ) and m^n. 
If rank (A) = n, where 

I dx\ dx„ 
dyi ' ' ' dyi 



(13) 



A 



\ 



dxi 
dy m 



dx„ 
dy m 



then V/(x) = if and only if Vit(y) = 0. 

Let us now consider the 3-body case, and for this purpose we will 
use the notation rj = (^1,^2) for i = 1,2,3. From Lemma H we have 
that if rank (A) = 3, where 



A 



dri2 
dqu 

dri2 
dqi2 
dri2 

9<?21 

dri2 
9<?22 

dri2 

9Q31 

dri2 

9<?32 



dri3 
dqu 

9ri3 
dqi2 
§ri3 
dq2l 
dri3 
dq22 
drn 
dqu 
dris 
dq-32 



dr 2 3 
dqu 

dr 2 3 
dqi2 

dr 23 
9<J2i 

dr 2 3 
9(J22 

dr23 

9<?3i 
dr 2 3 
dqz2 



\ 



I 



qu -921 


911-931 





ri2 


m 


912-922 


912-932 





ri2 


ri3 


911-921 





921—931 


ri2 


r-23 


912—922 





922—932 




f23 





911-931 


921-931 


ri3 


r23 





912—932 


922—932 


ri3 


»*23 
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then V[/(ri,r 2 ,r 3 ) = if and only if Vf7(?r 2 , ?T 3 , r 23 ) = 0. 

Some straightforward computations show that the rank (A) = 3 if 
and only if 

qu qu 1 
det | q 2 i q-n 1 | + . 

<?31 932 1 

This determinant is twice the oriented area of the triangle formed 
by the 3 particles. In short, if r l5 r 2 and r 3 are not collinear, then 
W (r 12 , r 13 , r 23 ) = if and only if W (n, r 2 , r 3 ) = 0. 

Using Lemma H] in order to find the planar central configurations, we 
first need to solve the equation 

(14) W = aVI 

in terms of the mutual distances r^, taking into account the fact that 
the moment of inertia, J, can be written in terms of the mutual dis- 
tances as / = (1/m) Xir=i m i m j r iji where ifi is the total mass. So, for 
fixed % and j, we have 

rriimj miirij a 

— b — rr-r- = 2 — m,'m,-r. 



r 2 " 6+1 tj- 

' ij r ij " L 

Multiplying by r\+\ we obtain 

/(r fi ) := 2ar b + 2 + rhr^ 1 + mb = 0. 

Regarding the above equation as a polynomial in the variable r^, 
since a < 0, /(0) = bm > and the coefficients polynomial have 
just one change of sign, we can verify easily that the function / has 
exactly one positive root. Observe that the function / only depends 
on the total mass rh, and therefore the respective solution for /(r^) is 
the same for all mutual distances. We have thus proved the following 
result. 

Theorem 2. In the Manev-type three body problem, for any values 
of the masses, there are exactly two equilateral central configurations, 
which correspond to the two possible orientations of a triangle in a 
plane. 
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6. A Framework for the Study of Collisions 

We will further study the dynamics at an near total collision for 
Manev-type n-body problems. A convenient framework for this pur- 
pose is given by the so-called McGehee coordinates [8], 

p = (r'Mr) 1 / 2 
s = p _1 r 

(15) v = p b /\p t s) 

u = / /2 (p-( P < s)Ms), 

where M = diag(m 1; m 1; m 2 , m 2 , . . . , m n , m n ), r = (ri,...,r n ), and 
p = (pi, . . . , p n ). After a reparametrization of the time variable, 

(16) dr = r~ l - b/2 dt, 

the equations of motion (jSJ) become 

p' =pv 
b 



v 1 =- v 2 + u*M _1 u - p 6-1 W(s) - bV(s) 



(17) 



u' 



=M _1 u 

'b 

„6-l 



U - (u*Af _1 u)Ms+ 



p"" i [iy(s)Ms + W(s)] + 6y(s)Ms + W(s). 

Here the prime denotes differentiation with respect to the new (ficti- 
tious) time variable r, and the old notation is maintained for the new 
dependent variables, which are now functions of r. Furthermore, the 
new variables fulfill the constraints s*Ms = 1 and u's = 0. 

In these coordinates the energy integral ([6]) turns into the relation 

(18) ^(u f M _1 u + v 2 ) - p^Wis) - V(s) = hp b . 
We define the total collision manifold as 

(19) c = {{p,s,v,u)\ p = 0, u*M -1 u + i; a -2V(s) = 0}. 

Notice that p 1 = if p = 0, so C (which is an analytic submanifold of 
codimension 1 in the boundary of the phase space) is invariant under 
the flow of the system (fT7|) . By continuity of the solutions respect to 
initial conditions, the flow on C provides important information about 
the orbits close to triple collision (see [8] for more details). The total 
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collision manifold can also be regarded as an invariant boundary pasted 
onto each energy surface: 

(20) E h = {(p, v, s, u^u'M-'u + v 2 ) - p b - l W{s) - V(s) = hp b }. 

These concepts are ideal for understanding the qualitative behaviour 
of total- and near-total-collision solutions. 

7. Collision and Near-Collision Dynamics 

In this section we will study the dynamics of total- and near-total- 
collision orbits of the Manev-type n-body problem. An important role 
in this study is played by central configurations and by the solutions 
that can be derived from them. 

In the planar Newtonian n-body problem, a rigid rotation of a central 
configuration is called a relative equilibrium; in rotating coordinates, 
relative equilibria are fixed points. A non-rotating homothetic orbit of 
a central configuration is called a homothety. The composition of a 
relative equilibrium and a homothety is called a homographic solution. 

In the Manev-type three-body problem, since the potential only de- 
pends of the bodies' mutual distances, the central configurations are 
invariant under rotations, so any central configuration determines a 
particular periodic orbit, which in a rotating frame is a fixed point. 
So in the Manev-type three body problem, any central configuration 
corresponds to a relative equilibrium. 

In the Newtonian case, any central configuration also corresponds to 
a homothetic orbit. But is this valid for Manev-type potentials too? 
As we will further prove (see Theorem [5] and Section 8), this property 
is not satisfied in general. To show this, and to determine under what 
circumstances homothetic solutions still exist, we will prove several 
preliminary results. 

Notice that the flow on C is given by the equations 

b 

v > = - v 2 + u*M u - &V(s) 

(21) s' = Af -1 u 

u' = (^ - l)uv - (u t M" 1 u)Ms + bV(s)Ms + W(s). 

The equilibrium points of system (T2T]) are given byu = 0, v = ±^/2V(s), 
where s must be a critical point for the function V(s) restricted to the 
unit sphere corresponding to the mass matrix M. The masses are in- 
volved because the equation bV(s)Ms + VV(s) = must be satisfied. 
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But these are the critical points of the function V, which is the re- 
striction of the homogeneous potential V to the unit sphere given by 
the mass matrix M. Such critical points correspond to the central 
configurations of the homogeneous potential V. 

Proposition 3. For any value ofb>2, the flow on the total collision 
manifold C is gradient-like with respect to the coordinate —v (i.e. the 
flow increases with respect to —v along non- equilibrium solutions). 

Proof. The energy relation (fl8|) . restricted to C, takes the form 

u*M _1 u + v 2 - 2V(s) = 0. 

Using the above expression and (12"T|) . we get that 

(22) v' = (l- ^)u*M _1 u 

on C. If u 7^ then v' < is increasing with respect to —v. On 
the other hand if u = then u' = bV(s)Ms + W(s), i.e. u' = 
only if s is a critical point of V. Consequently — v is strictly increasing 
along nonequilibrium solutions, which means that the vector field is 
gradient-like with respect to —v. □ 

Denote by ind(s ) the index of the critical point s , i.e. the number 
of eigenvalues of D 2 V(sq) with negative real part. Then we can prove 
the following result. 

Theorem 3. Let s be a nondegenerate central configuration of the 
planar n-body problem with potential V and b > 2. Then the dimensions 
of W u (sq ) and W s (sq) are the same and equal to 2n — 2 — ind(s ) in 
Eh- The dimensions ofW s (sQ ) and ^"(sq ) are the same and equal to 
2n — 4 + md(so) in Eh- The dimension of Eh is 4n — 5. 

Proof. Let s be a central configuration, v = ±a/21 / (s ), and u = 0, 
then the equation of motion restricted to Eh are 



p =pv 

1 - 

2 



v ' = [\- b \ u'M^u + (b- l)p h - l W(s) + bhp h 



(23) s' =M- X u 



u' = ( - - 1 ) uv - (u t M- 1 u)Ms+ 



2 

p b - 1 [W(s)Ms + W(s)] + bV(s)Ms + VV r (s). 
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Taking into account the centre of mass and linear momentum inte- 
grals as well as the restrictions s t Ms = 1 and u*s = of the McGehee 
coordinates, the above system has dimension An — 4. 

Linearizing the system, the eigenvalues for 6 > 2 are given by the 
matrix equation 

/ v ... ... \ 

* ... ... * 



(24) 



A (| - l)vl 2n . 



— 



4n— 4 



o 



4n-4; 



/ 



\ o o 

where In is the N x N identity matrix, On is the N x N zero matrix, 
A denotes the Hessian matrix of V (i.e. the potential restricted to the 
sphere of constant moment of inertia) and * denotes an element without 
importance in the computation of the eigenvalues. 

It is clear that the first two eigenvalues are v ^ (since V(so) ^ 0) 
and 0. To obtain the remaining eigenvalues of equation (1241) . suppose 
z is a (2n — 3)-vector satisfying 

Az = A,z 



(25) 

for i — 1, . 
A. Then 

02n-3 

A 



2n — 3,i fixed, where Ai, 



j A2n-3 



are the eigenvalues of 



^2n-3 \ ( z 

(6/2 - l)vI 2n -3 )\^z 
Consequently \i is a root of equation (T54 



{Ai + (6/2 - 



if 



1? - (6/2 - !>// - Ai 



which gives 



!' 2 = ^{(6 - 2)v ± ^(2 - 6) V + 16Ai} 



2n-3. 



for i = 1, 

Then if t> = a/^(s ), the differential matrix of the vector field re- 
stricted to Eh has 2n— 2— ind(s ) eigenvalues with positive real part and 
2n — 4 + ind(so) with negative real part. The values of the dimensions 
are switched if v = — a/V(s ). □ 

Theorem 4. Le£ So 6e a central configuration of the collinear n-body 
problem with potential V and 6 > 2. Then the dimensions o/ ^"(sj) 
and W s (sq) are the same and equal to n — 1 in Eh- The dimensions 



of W s (s+) and W u (s 
dimension of Eh is 2n 



are the same and equal to n 
3. 



2 in E h . The 
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Proof. The proof is similar to the one of the previous theorem. Let Sq 
be a central configuration, v = ±y/2V(so). The equations of motion 
restricted to Eh are given by equation (|23|) . with the obvious modifica- 
tions. 

Linearizing the system, the eigenvalues in the case b > 2 are given 
by the following matrix equation 

/ v ... ... \ 

* ... ... * 



(26) 



O n -2 

A 



\ o o 



In-2 

- lW n _ 2 



2n-2 



o 



2n-2i 



I 



where In and On are defined as before. Again A is the Hessian matrix 
of V and * denotes an element without importance in the computation 
of the eigenvalues. The first two eigenvalues are v ^ (since V(sq) ^ 0) 



and 0. If Ai, 



_ 2 be the eigenvalues of A, then 



/4' 2 = - A {(b - 2)v ± y/(2 - 6) V + 16A,} 

for i — 1, . . . , n — 2. Note that, in this case, Lemma 4 implies that A is 
positive definite, and thus the eigenvalues Ai, . . . , A n _ 2 are all positive. 
Consequently, for v > 0, fij is negative and fif positive. However, 
for v < 0, fij is positive, whereas /if is negative. This concludes the 
proof. □ 

We will further state and prove a result that clarifies under what 
circumstances homothetic solutions exist. 

Theorem 5. A solution of the Manev-type n-body problem is homo- 
thetic if and only if the particles form, at all times, a simultaneous 
central configuration for the potentials V and W . 

Proof. Assume that the solution is homothetic, then s = s , where s 
is a constant. Therefore s' = and, from the second of equations (|2"T|) . 
u = 0. Thus the homothetic orbits are confined to the invariant plane 

(27) P = {(p,s,i;,u)|s = So,u = 0}. 

So u' = implies that p b - 1 [W(s )Ms + VW(s )] + bV(s )Ms + 
Vy(s ) = 0. If p is not constant then there are Pi ^ and p 2 ^ with 
Pi 7^ P2 such that 

p^ 1 [H/(s )Ms + W(s )] = -[bV(s )Ms + VV(s )] 

p^ 1 [W(s )Ms + W(s )] = -[&V(s Q )Ms + W(s )]. 
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This means that [M/(s )Ms +W(s )] = and [W(so)Afs +VW(s )] = 
0, i.e. that s is a simultaneous central confiiguration for the potentials 

V and W. If p is constant, p' = and either p = or v = 0. The first 
case is trivial, whereas in the latter case — p 6_1 VT(s) — bV(s) = 0. But 
this is impossible since V > and > 0. 

If s = So is, at all times, a simultaneous central configuration for V 
and W, then the solution is obviously homothetic. □ 

The next result proves the existence and uniqueness of heteroclinic 
homothetic solutions for b > 1. 

Theorem 6. Let s be a simultaneous central configuration for the 
potentials V and W . Then, if b > 1, every energy surface of nega- 
tive constant (h < 0), contains a unique homothetic solution defined 
on (— oo, oo), satisfying s = s for all times and such that p(r) — > 
when t — ► ±oo. In other words the solution begins and ends in a total 
collapse, maintaining for all times the same central configuration. 

Proof. Since So is a simultaneous central configuration for the potentials 

V and W, we have that [bV(s )Ms + W(s )] = and [W(s )Ms + 
VW(s )] = 0. Consequently, the set 



is invariant for the equations of motion. Restricting these equations to 
V, we get 



V = {(p, s,t>, u)|s = s , u = 0} 



P 



= pv 



(29) 



' = ^2_ p 6-l^ (so) _ 6F(so); 



V' 



while the energy relation becomes 




W(so) - V(so) = hp b . 



Equations fl29|) become 



P 



pv 



v' = (b 



l)p b - l W(s ) + bp' 



< b h. 



This leads to 



dv 1 
dp v 



Kb 



l)p b - 2 W(s ) + bp b - l h] 



which yields 



(30) 



2 



p b - l W(s ) +p b h + K. 



16 



F. Diacu, E. Perez-Chavela and M. Santoprete 



where, if b > 1, we choose K = V^(so). If b > 1 and h > 0, then 
\v\ > ±a/2V(s ) and the homothetic orbits are not heteroclinic. If 
h < 0, there is a unique curve connecting the points (a/2V(so), 0) and 
( — a/2V(s ), 0) on the plane (v , p). These facts prove the theorem. □ 

The following result shows that the above property is also true for 
the equilateral central configurations. 

Corollary 1. Let s be an equilateral central configuration for the po- 
tential U. Then, if b > 1, every energy surface of negative constant 
(h < 0) contains a unique heteroclinic homothetic solution. 

Proof. Clearly s is a simultaneous central configuration for the poten- 
tials V and W. The proof follows from Theorem O □ 

Two submanifolds £\ and £2 of a submanifold £ are said to be trans- 
verse at a point x if one of the following situation arises: 

(1) £ 1 n£ 2 = $; 

(2) x G £\ fl £2 and T x £\ + T x £ 2 = T x £, where T x £ denotes the 
tangent space to £ at the point x. 

We can now prove the following result: 

Theorem 7. In the planar Manev-type n-body problem with b > 2, a 
necessary condition for having a transversal homothetic solution 7h(s ) 
in Eh with h < is that V be a non- degenerate minimum at the point 
So associated with the homothetic solution. 

Proof. Let 7h(so) be a transversal homothetic solution in Eh with h < 0. 
Then by Theorem 3 both W^(s + ) and ^(s^) are (2n-2-ind(s ))- 
dimensional and Eh is (4n — 5)-dimensional. Since 7h(so) £ ^"(so) D 
W s (sq) and 7h(sn) is transversal we have that 

dim^ < dimiy u (s ) +dimVy s (s ) - 1. 

That is An — 5 < An — 5 — 2 ind(s ). Therefore ind(s ) = and the 
function V has a nondegenerate minimum at s . □ 

8. Simultaneous Configurations and Relative Equilibria 

In this closing section, we will show that, for most choices of the 
masses in the quasihomogeneous 3-body problem, the collinear central 
configurations of the potential U are not simultaneous relative equilib- 
ria for V and W. 

Theorem 8. Let S3 be the set of masses (777.1,7712,7713) G for which 
the collinear configurations are simultaneous central configurations for 
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the potentials V and W . Then the set E 3 is nonempty and nowhere 
dense in . 

Proof. Assume the configuration sy (mi, 777,2, TO3) is a collinear central 
configuration for V and S\y(mi,m 2 ,m 3 ) is a collinear central configu- 
ration for W. In [7j, Euler found a complicated formula that expresses 
the ratio of the distances between the masses for any rectilinear central 
configuration in the Newtonian case. Euler's formula can be directly 
extended to any homogeneous potential. Moreover, the fact that Eu- 
ler's expression is an analytic function of the masses remains true in the 
homogeneous case. Therefore both Sy and S\y are analytic functions of 
mi, m2 and 1713, as long as the masses are positive. Consequently the 
function z = Sy — S\y is also an analytic function of the masses. 

For the function V, Euler's formula depends on a, whereas for W 
it depends on b. So in general Sy 7^ therefore for every a and b 
with a 7^ b there are values of the masses for which z 7^ 0. Since z is a 
nonzero analytic function, its zeroes form a nowhere dense set. 

The nonemptiness of the set of simultaneous central configurations 
follows from noticing that if mi = m 3 and the mass m,2 is located 
halfway between the other two, then the three masses form a simulta- 
neous central configuration for V and W. 

□ 

A consequence of Theorems |5] and [8] is that, for most values of the 
masses, there are no rectilinear homothetic orbits. More precisely: 

Corollary 2. If (mi, 777.2, ^3) £ \ ^3, then there are no rectilinear 
homothetic orbits. 

This shows that the rectilinear homothetic orbits are characteristic 
to homogeneous potentials, but they prove unlikely in the quasihomo- 
geneous case. 
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